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ABSTRACT 

We investigate the cross spectrum of XTE J1701-462 in various types of 
Neutron Star-Low Mass X-ray Binary subclasses, during its 2006-2007 outburst. 
We analyze the relation between the time lags and temporal variabilities. We 
find that the hard time lags accompany with the horizontal branch oscillations 
(HBOs) and the soft time lags dominate the noise in the low frequency range 0.1- 
10 Hz on HB. In the Cyg-like phase, the time lags decrease on the middle normal 
branch (NB) from HB/NB vertex to NB/FB vertex, whereas the time lags are 
roughly invariant in the Sco-like source. We discuss that the Compton upscat- 
tering by the corona introduces the soft lag in low frequency noise. We interpret 
that the variation of the Comptonization component from the disk emission lead 
to the HBOs' time lags evolution along the Z tracks. We also report the rms 
amplitude spectrum and phase lag spectrum for the NBO. The ~ 160° phase 
lag is found. We present that the rms amplitude of both the Cyg-like and the 
Sco-like NBO linearly increase with the photon energy in low energy bands, and 
it will drop in the highest energy band. 

Subject headings: binaries: general — stars : individual (XTE J1701-462) — stars: 
neutron - X-rays: binaries — X-rays: individual (XTE J1701-462) — X-rays: stars 
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Introduction 



The Low mass X-ray binaries (LMXBs) are binary syste ms which conta i n a N eutron 



van der Klis 



fl2006f ) for 



Star accreting material from a low mass companion star (see 
review). Based on the correlated X-ray spectral and timing properties, they can be divided 
into two subclasses, the Z sources with near Eddington luminosity (~ 0.5 — 1 L e dd) and the 
atoll sources with much lower luminosity (~ 0.0001 — 0.5 L e dd)- The Color-Color Diagram 
(CD) and Hardness-Intensity Diagram (HID) of Z source have three branches, from top to 
bottom, which are so-called the horizontal branch (HB), the normal branch (NB) and the 



flare Branch (FB) ( IHasinger &: van der Klis 



19891 ). respectively. However, the CDs/HIDs 



of atoll sources are obviously different from those of the Z sources. A CD/HID of atoll 
sources also have three branches, the extreme island state, the island state, and the banana 
state. Usually, a complete track of the Z sources last from hours to days, whereas atoll 
sources last from weeks to months. 

In terms of the orientation of branches, the Z sources are classified into the Cyg-like Z 
sources (Cyg X-2, GX 5-1, and GX 340+0) with a horizontal HB ("Z"-shaped tracks) and 
the Sco-like Z sources (Sco X-l, GX 17+2, and GX 349+2) with a vertical HB ( V-shaped 
tracks). The Fast Fourier Transform (FFT) algorithm was applied in order to obtain the 
power density spectra (PDS) of LMXBs. In some Z sources (Cyg X-2,GX 5-1, GX 340+0, 
GX 17+2, and so on), the HB Quasi-periodic oscillation (QPO) signals, which vary with 
frequencies in range from 13 Hz to 55 Hz, are observed in the power spectra. Almost all 
Z sources have 5-7 Hz QPO on NB and few of them have ~ 20 Hz QPO on FB. The 



radiation hydrodynamic model ( iLamb 



1989 



Fortner et al. 



1989 



Miller fc Lamb 



19921 ) 



was proposed to explain the ~ 6 Hz NBOs. In this model, the NBOs of Z sources arise 
from the optical depth oscillation of in-flowing medium near the neutron star scatting by 
the outgoing X-ray flux. This model also explained a minimum in the NBO rms amplitude 
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spectrum when the phase lag spectrum showed 180° phase jump. In C yg X-2, a 150° phas e 
jump was observed aroun d 6 keV, which was defined as pivot energy (IMitsuda fc Dotani 



1989 



Dieters et al. 



2000). Moreover, the same phase jump was observed around 3.5 keV 
in GX 5-1. Both of these two Cyg-like sources displayed a minimum around pivot energy 
in the NBO rms amplitude spectrum. However, the Sco X-l which belongs to Sco-like Z 
sources did not emerge a phase jump on the NB. The near-180° phase jumps observed in 
Cyg X-2 and GX 5-1 but not in Sco X-l are consist with the classification of Z sources. 



Two kinds of time lags were observed in LMXBs. Firstly, iLei et al. I ( 120081 ) applied 
the cross correlation function to Cyg X-2 during its low-intensity states, and detected 
both anticorrelation and positive correlation with hectosecond soft and hard time lags. 
he anticorrelation and time lags were also observed in HB and upper NB of GX 5-1 



Sriram et al. 



van der Klis et al. 



2012) whi ch implied a truncated accretion disk in LMXBs. Secondly, 



(119871 ) used the Fourier cross spectrum technique to study the Cyg X-2 
and GX 5-1. They found tens of milliseconds hard time lags, i.e., the har d photons lag 



the so ft photons in HBOs and NBOs, and soft lags in low-frequency noise. 



(119991 ) discovered millisecond soft lags in the atoll source 4U 1636-536. 



Kaaret et al. 



Qu et al. 



(2001) 



found that Cir X-l presents a hard lag on the HB and a soft lag on the NB and FB. In 



order to exp 



ain both soft lags an d hard lags in LMXBs, several mech anisms, inc 



shot model (jVaughan et al. 



two-layer corona mo del (jNobili et al. 



1994 ), a uniform Comptonizati on model (ICui et al 



2001 



(jUttlev et al. 



Qu et al. 



uding a 



1997al ). a 



20011 ) and a disk propagation model 



201 lh . had been proposed. 



The XTE J1701-462 is a unique LMXB sourc e observed 



Explorer (RXTE) during its outburst in 2006-2007 (IRemillard 



jy Ro ssi X-ray Timing 



20061 ). At high luminosity, 



the XTE J1701-462 switched between a Cyg-like Z source and a Sco-like Z source and at 



low luminosity the transition from Z source to atoll source was observed (IHoman et al. 
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20071 : iLin et al. 1120091 : iHoman et al. 1 12010( 1 . During these states, the shape of CDs/HIDs, 
the X-ray timing behavior and the X-ray spectra also changed. It provides us a great 
opportunity to understand the NBO behaviour, the time lag spectrum and the rms 
amplitude spectrum, between Cyg-like phase and Sco-like phase. We also utilize the cross 
spectrum to study the time lags evolution along with branches from Z source to atoll source. 

In Section 2, we describe the observations and analysis techniques and our results are 
reported in Section 3. In Section 4 and 5, we present the discussion and conclusions. 



2. Observations 

We analysed all 866 observations of the LMXB XTE J1701-462 from its 2006-2007 
outburst collected with the Proportional Counter Array (PCA) instrument on board of 
RXTE. During the outburst, the source showed X-ray spectral and timing behaviours of 
the Z source and the atoll source. 

For all observations, we employed the following standard data selection criteria, a 
source elevation of > 10°, a pointing offset of < 0.01° and a SAA exclusion time of 30 
minutes. We utilized the Heasoft 6.12 to extract background-subtracted light curves with 
a 16 s time resolution from the "Standard 2" mode data and no dead-time correction was 
applied. Only data from PCU 2 were used. In order to obtain the CDs/HIDs from the 
light curves, we define the soft color as 4.5-7.4 keV/2. 9-4.1 keV (channels 10-17/6-9) count 
rates ratio, while the hard color as 10.2-18.1 keV/7.8-9.8 keV (channels 24-43/18-23) count 
rates ratio, then defined the intensity as count rates covering the energy range 2.9-18.1 keV 
(channels 6-43). Our data selection is listed in Table Q] and the CDs/HIDs are displayed in 
Fig. [TJ The interval I and II belong to the Cyg-like source, Interval III and IV belong to the 
Sco-like source. 
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We used the single-bit and science event model data with 32 s segments and 2 -9 s time 
bin to created power density spectra (PDS) in frequency range 1/32 Hz-128 Hz. We studied 
the HBOs in energy range ~ 2 — 60 keV and the NBOs in energy range ~ 7 — 60 keV, 



becau se the NBOs were more apparent in ~ 7 — 60 keV than in ~ 2 — 60 keV (IHoman et al. 



20071 ). No background-subtracted an d no dead time correctio n were applied. We employed 
the power spectral normalization of iMivamoto et al. 



power spectral using a multi-Lorentzian fun ction (Nowak 



Inter val I- IV. All parameters were set free (jPottschmidt et al. 



(119911). And then, we fitted the 



2000 



Belloni et al. 



2003 



200 21) for th e 



van Straaten et al. 



20051 ). We obtained the centroid frequency v c , the full width half-maximum (Az/, FWHM) 
and the fractional root-mean-squared (rms) of HBO. We plotted the power spectra in the 
power times frequency representation (vP v , in units of rms 2 ), where each power spectral 
density P v is multiplied by its Fourier frequency v. 

We also extracted five energy bands, i.e., 2.06 - 4.49 keV (channels 0-10, CHI), 4.9 - 7 
keV (channels 11-16, CH2), 7.4 -9.4 keV (channels 17-22, CH3), 9.8 - 14.8 keV (channels 
23-35, C H4), 15-65 keV (ch annels 36-149, CH5), to compute cross spectra. The cross 



spectrum (INowak et al. 



1999J) is defined as C(j) = S(j)*H(j), where S(j) and H(j) are 



the measured complex Fourier coefficients at a given frequency uj for soft and hard energy 
bands, respectively. The phase lag between two energy bands light curves is the Fourier 
phase 4>(j) = arg[C(j)]. The time lag is constructed from <p(j) by dividing by 2iTVj, i.e., 



T[Va 



<fi(j) /2iTVj. A positive time lag indicates that the hard photons lag the soft photons. 



The cross spectra of XTE J1701-462 show that the phase lag above 100 Hz are consistent 
with 7r. This is dead time eff ect which should be corrected by subtracting a cross vector, 



averaged over 100 to 128 Hz ( Ivan der Klis et al. 



1987 



Vaughan et al. 



QPO s and kHz QPOs were observed in XTE J1701-462 (IHoman et al. 



2010 



1999) Hectohertz 



Sanna et al. 



20101 ). the subtracted cross vector might be combined with the white noise and signal. 



However, the signal above 100 Hz can be dismissed compared with the noise amplitude. We 
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computed the time lags of CH2 - CH5 relative to CHI of NBO and investigated the time 
lags evolution along with HB between 9.8 - 14.8 keV and 2 - 4.5 keV. 



To investigate the correlation between the data p osition in the HID and the temporal 



properties, we introduce a rank number S z defined by 



Hasinger et al. 



fll990h . In Fig. El we 



set the S z of box 12 and 24 as 1 and 2, respectively. The rank numbers of oth er boxes are 



decided by spline interpolation (IDieters fc van der Klis 



2000 



Lin et al. 



20121 ) 



3. Results 

3.1. Time lags evolution along HB in the Cyg-like phase 

We calculate the PDSs and average cross spectra for XTE J1701-462 in the Cyg-like Z 
phase (Fig. [3] and Fig. HJ) and in the Sco-like Z phase (Fig. |5]) . Positive lags indicate that the 
hard photons (9.8-14.8 keV) lag the soft photons (2-4.5 keV). In Fig. [31 the cross spectra 
of Interval I show hard lags (red '+') during QPOs' frequency range on HB. The red noises 
below 10 Hz display soft time lags (black 'o') on HB. The time lags of NB are dominated 
by hard lags. In Fig. HJ the HBOs and cross spectra of Interval II are present. The time 
lags on HB of Interval II and Interval I have similar behaviours. In Fig. [5j the hard lags are 
accompanied by QPOs and the soft lags are observed less than 10 Hz on HB. Meanwhile, 
we plot the phase lag versus frequency in Fig. [3H5J which clearly show the intrinsic phase 
lag couple with QPO appearing. 

In order to investigate the correlation between S z and temporal variabilities, we average 
the time lags of QPOs over the frequency between v c — Au/2 and u c + Au/2. In Fig. |6l we 
display the QPO centroid frequencies, the time lags and rms of the HBOs as a function of 
S z in Interval I. The centroid frequencies of HBO increase from 12 Hz to 53 Hz, but the 
time lags and rms of HBO decrease when S z is below 0.5, then remain constant when S z 
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is larger than 0.5. As a contrast, we provide the same temporal variabilities of Interval II 
in Fig. [7J We notice that the centroid frequencies and the rms of HBOs are similar with 
interval I. However, The time lags rapidly decrease when S z less than 0.5, then the time 
lags slowly increase with S z . All time lags are hard lags with several milliseconds. 

For different HIDs, the intensity of the HB/NB vertex and the NB/FB vertex were 
considered as reference value for the accretion rate fa. In order to study the time lag - 
accretion rate connection, we also calculate the time lag between 7.4 - 9.4 keV and 2 - 4.5 
keV in the HB/NB vertex region of Interval I- IV. In Fig. [HJ the QPO centroid frequency 



positively correlates with the time lag. Both of them are 



ollowe d by source intensity 



decreasing, that is the accretion rate decreasing (jLin et al. 1 120091 ) . We should notice that 
the source intensity difference between Interval I and Interval II are small (upper panel in 
Fig. [1]) and the time lags are extremely close too. 



3.2. Time lags evolution of the Z phase on NB 

In Fig. [9j we show the variation of NBO in Interval I. The centroid frequencies of NBO 
are around 6 Hz which consist with other Z sources. In order to study the time lags along 
NB, we average the time lags with a fixed frequency range 6±2 Hz, between 7.4 - 9.4 keV 
and 2 - 4.5 keV. 

In Fig. [TUJ we display the time lags evolution from HB/NB vertex to NB/FB vertex of 
Interval I. The time lags are nearly constant on upper NB (1 < S z < 1.2) and lower NB 
(1.7 < S z < 2), and reduce from ~ 45 ms to ~ —70 ms on middle NB (1.2 < S z < 1.7). 

In Fig. [TTj the time lags as a function of S z on NB and FB of Interval III are reported. 
Between HB/NB vertex and NB/FB vertex, the time lags averaged over 4-8 Hz are nearly 
constant around 50 ms. 
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3.3. The time lags spectrum and rms amplitude spectrum of NBOs in the 

Cyg-like phase and the Sco-like phase 

We extracted the the Cyg-like phase and the Sco-like phase NBOs with five energy 
bands light curves, 2.06 - 4.49 keV (channels 0-10, CHI), 4.9 - 7 keV (channels 11-16, 
CH2), 7.4 -9.4 keV (channels 17-22, CH3), 9.8 - 14.8 keV (channels 23-35, CH4), 15 - 65 
keV (channels 36-149, CH5), respectively. The NBOs' time lags of CH2 - CH5 against 
CHI are computed. For each energy band, the rms of NBO is obtained by fitting the 
relative PDS with a fixed centroid frequency and FWHM. We defined the time lags as a 
function of energy "time lag (or phase lag) spectrum" and the rms as a function of energy 
" rms amplitude spectrum" . 

In Fig. [T21 we notice that the maximum time/phase lag on the lower NBO is 55±7 ms 
or 114 ± 15° in Interval I, respectively. On the lower NB of Interval III , the time/phase lag 
is 49±6 ms or 106 ± 13° (Fig. [T3"|) . On the lower NB of Interval IV, the time/phase lag are 
slightly larger than Interval I and III, i.e., 74±16 ms or 160 ± 13° (Fig. [T4"j) . The near-160° 
phase jumps occur around 10 keV in both the Cyg-like phase and the Sco-like phase. 

Regardless of CH5, the rms amplitude of NBOs linearly correlate with energy in the 
range 2-15 keV (i.e., the signal strength of QPOs are proportional to the photon energy.) 
of Interval I and III (right pa nels in Fig. [T2l and Fig. [TBI . The rms energy dependence 



was also observed in Sco X-l (jWang et al 



20121 ) . Because of low count rates and low 



signal-to-noise in the highest energy band, the rms of NBOs have large error bars. The 
drop of the rms amplitude in the highest band was appeared. A pivot energy in 4.9-7 keV 
appears in the rms amplitude spectrum of Interval IV (right panel in Fig. [T41) but with a 
weak tendency. 
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4. Discussion 



The LMXB source XTE J1701-462 transited from a Z source into an atoll source during 
its outburst. We used the FFT to compute the PDS of XTE J1701-462. For the Cyg-like 
observations (Interval I & II), the centroid frequencies of HBO increase from 12 Hz to 53 
Hz, while the amplitude of HBO becomes weaker. The cross spectra of HBO were computed 
as well. In Fig. [3] and Fig. HJ we found that the time lags of low frequency noise less than 
10 Hz are dominated by soft lags. Moreover, the soft lag increases towards low frequency. 
When QPOs appeared, the XTE J1701-462 showed hard lags, i.e., high energy photons 
lag soft photons. For the NB observations of the Cyg-like phases and the Sco-like phases, 



the hard lags w ere also observed during NBOs' frequency 



other Z sources (Ivan der Klis et al. 



1987 



Vaughan et al. 



rang e . These a 



1994! ; 



Qu et al- 



so appeared in 



2001 



200J). We 



also study the time lags evolution along NB in the Z sources. In Interval I, the time lags 
averaged between 4 Hz and 8 Hz decrease on middle NB as a function of rank number S z 
(Fig. [TO]) . This could be the centroid frequency and rms amplitude variation of the NBO 
(Fig. [9]). In interval III, the time lags of NBO nearly remained 50 ms (Fig. [TTlh For the 
HB/NB vertex regions in the Z sources, the time lag increased with the source intensity 
decreased, which is the trend of accretion rate decreased. 

In order to explain the observed time lag in compact star system(e.g., BH-, NS-LMXBs 



AGNs, and so on), many models are proposed, including the shot model (Vaughan et al. 



19941 ). the propagation o 



Comptonization models (ICui et al 



accretion fluctua t ions model (lUttlev et al. 



1997b 



Kazanas et al. 



1997) 



20111 ). and the 



mathematically explained the observed power spectra in LMX Bs (lAlpar fc Shaham 



The shot mode 



128a), 



the soft lag in t he low frequency noise (IShibazaki et al. 1119871 ) and the hard lag in the HBO 



simultaneously (IVaughan et al. 



19941 ) . It regarded the HBO as the signature of intensity 



modulation and the low frequency noise as the power spectral signature of shots, when the 
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accretion matter fall onto the Neutron Star surface. The observed soft lag of low frequency 
noise could be interpreted if the shot envelopes in the soft energy band lag those in the 
hard energy band. Meanwhile, the hard lag in the HBO may originate from a particular 
delaying of the QPO production in high energy band without distinctly interfering by low 
frequency noise. 

The time lag in the hard state of BH-LMX Bs which displays ha r d lag in low frequen cy 



noise shows the exactly 



2003 



Uttley et al. 



o pposite to our r esults (INowak et al. 



1999 



Pottschmidt et al. 



201lh . 



Nowak et al. 



(Il999h investigated the hard state time lag in the 



Cyg X-l and found that the energy dependence of hard lag in low frequency noise as well 
as the hard lag increase with decreasing low Fourier freq uency. It seem s cont radictory to 



Uttley et al. 



(1201 lh pointed out 



the uniform and compact Comptonization corona model, 
that the observed hard lag in GX 339-4 corresponds to the viscous propagation of mass 
accretion instabilities through the disk. The emission of BH-LMXBs contains two parts, 
the power law dominant in the high energy band and the disk black body dominant in the 
low energy band. When the mass accretion fluctuations happened, the variabilities of disk 
black body emission lead to the power law emission variations on time scales of seconds. 
If the low frequency noise comes from these two components variations, then the hard lag 
appears. For the soft lag, it may be responsible by the power law illuminating the black 
body on time scales of milliseconds, where the the power law emission regions locate about 
tens of Rq to the disk. 



How ever, XTE J1701-462 shows more 



LMXBs ( Lin et al. 



2009 



Ding et al. 



complicated emission mechanism than BH- 



20111 ). On the HB, the energy spectral of XTE 
Jl 701-462 is composed of the multi-color black body from the disk, the black body from 
the NS surface and the cutoff power law from the corona above the disk, which dominate 
the emission in several keV, tens of keV and the hard tail respectively. On the NB and FB, 
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the cutoff power law component becomes too faint to be detected. We indica te that there 



exist three possible lag mecha nisms, the Compton upscattering by the corona (IKotov et al 



2001 



Arevalo fc Uttlev 



(lUttlev et al. 



20061 ). the propagation of the accretion rate instabilities in the disk 
201 if ) and the disk illuminating by the black body from the NS surface in 
XTE J1701-462. However, the soft lag introduced by the NS X-ray illuminating is likely the 
time scale of light travel from the NS to the disk, which is too low (~ 0.1 ms) compared to 
our calculations. The propagation of the accretion rate instabilities in the disk only produce 
hard lag in low frequency noise which is not conform to the soft lag in XTE J1701-462. 



The uniform Comptonization models (jPayne I Il980l ; iKazanas et al. 1 119971 ) can not 
interpret the energy dependence of hard time lag and soft time lags in LMXBs. In order to 
explain soft lags and time lags evolution of QPOs, a Com ptonization model with two-layer 



corona was proposed (jNobili et al. 



2001 



Qu et al. 



20011 ) . When photons pass through the 



two-layer corona, the photons encounter the inverse Compton scattering by hot electrons 
and gain energy. If the inner layer is hot and optical depth is < 1, the hard photons are 
scattering more than the soft photons and the hard time lags observed. If the inner layer 
is optical thick, the photons are efficiently upscattered. When the harder photons pass 
through the out part corona with lower temperature, they will lose energy and the soft lags 
appear. 



Qu et al. 



(j200lf ) explained that the accretion rate changing along branches leads to the 



Cir X-l time lag evolution. From multiwavelength campaigns of Cyg X-2, the monotonical ly 

~3 



increasing in rh from the HB, via the NB, to the FB was reported (lHasinger et al 



Vrtilek et al. 



199C 



1990|). If the accretion rate of XTE J1701-462 arise from the HB to the NB, 
the hot electrons in the inner layer corona are cooled by photons from the disk. This process 
will shrink the i nner lay corona and the observed HBO time lags bec ome smaller. Lin 



and co-workers (ILin et al 



2009 



Homan et al. 



2010 



Lin et al. 



20121 ) concluded that the 



13 



accretion rate did not vary significantly along the Z tracks but did along secular change, i.e. 
m decrease from the Z source to an atoll source. The time lag - QPO relation of the HB/NB 
vertex could be explained by the accretion rate decreasing (Fig. |H]). However, for the time 
lags evolution on the HB of XTE J1701-462, when the source intensity increased, the disk 
luminosity whi ch was mainly co ntributed by Comptonization ascended simultaneously 



(see Fig. 20 in iLin et al. 



(120091 ) ). When the luminosity of the disk component increases 
from the up turn to the HB/NB vertex, the temperature of the inner layer corona cool 
down through Compton emission and then lead to a smaller radius layer. Because the 
photons contributing to the HBO came from the disk, the time lags of HBO decrease as a 
function of S z was also observed. In a Z track, the time lag of HBO increasing towards high 
source intensity can be explained by the accr etion rate increasing and th e Comptonization 



component variations. 



Jackson et al. 



(120091 ) and 



Church et al. 



(120101 ) suggested an 



opposite accretion rate increasing direction, i.e., from the NB to the HB, which could not 
reasonably explain the time lags evolution of XTE J1701-462. 

The near- 180° phase lag and the dip in rms am plitude spectrum wer e observed in 



Dieters et al. 



(12000h provided two 



Cyg X-2 and GX 5-1, but did not show in Sco X-l. 
explanations for the NBOs behavior for Cyg X-2 and GX 5-1 and Sco X-l. Firstly, they 
may have different viewing angle. Secondly, with the radiation-hydrodynamic model for the 
QPOs, the QPO of GX 5-1 and Sco X-l are determined by optical depth oscillations while 
the QPOs of GX 5-1 are determined by luminosity variation, where the NS of GX 5-1 has a 
weaker magnetic field. XTE J1701-462 changes from a Cyg-like Z source, via a Sco-like Z 



source, to an atoll source in two years. The binary inc 



magn etic field are invariant on such a short timescale (IKuulers et al. 



ination and the strength of the NS 



1994 ; 



Psaltis et al. 



19951 ). The accretion rate dominates the whole processes. We found near-160° phase lag of 
NBOs both in Cyg-like phase and Sco-like phase which are consistent with Cyg X-2 and 
GX 5-1. However, no pivot energy in rms amplitude spectrum was observed in the Cyg-like 
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phase . This can not 



1989|). 



Wang et al. 



3e ex plained by the radiation-hydrodynamic model ( IFortner et al. 



( 120121 ) interpreted that the NBO and the rms energy dependence in 



Sco X-l originate from the transition layer. The similar behaviors of time lag spectrum 
and rms amplitude spectrum between the Cyg-like sources and the Sco-like sources indicate 
that the NBOs are irrelevant to the luminosity variation in XTE J1701-462. 



5. Conclusions 

We investigated the cross sprctral and the QPOs for the Cyg-like phase and the Sco-like 
phase of XTE J1701-462. The QPOs increase from 12 Hz to 53 Hz on the HB and remain 
~ 6 Hz on the NB. These consist with other Z sources. We also computed the time lags and 
HBOs for the HB/NB vertexes (Interval I- IV). The two-layer corona Comptonization model 
can explain the soft lag in low frequency noise and hard lag in HBO respectively. Both the 
accretion increasing and the Comptonization component variation can produce the HBOs' 
time lag evolution. For a NS-LMXB, we could not deduce its accretion rate increasing 
direction along branches in the HID only from the time lag variation. The spectral states 
should be analyzed simultaneously. This will be discussed in our future work. 
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Table 1. The data selection of the CDs and HIDs. 



Interval 


Begin of date 


Begin of Obs. 


End of date 


End of Obs. 


Source Type 




(DD/MM/YY) 




(DD/MM/YY) 






I 


21/01/06 


91106-01-04-00 


31/01/06 


91106-02-02-10 


Cyg-likc 


II 


02/02/06 


91106-02-02-14 


12/02/06 


91442-01-01-02 


Cyg-like 


III 


17/02/06 


91442-01-07-02 


25/02/06 


91442-01-03-03 


Sco-likc 


IV 


08/03/06 


92405-01-01-05 


15/03/06 


92405-01-02-05 


Sco-likc 
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Fig. 1. — CDs (upper panels) and HIDs (lower panels) for each of 4 intervals in table HJ 
Each dot represents 16 s of background-subtracted data from PCU 2. From left to right, we 
label Interval I- IV. 
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Fig. 2. — The HID of Interval I. We divide it into 26 boxes from up turn to FB. We set the 
rank number S z of boxes 12 and 24 to 1 and 2, respectively. The S z of other boxes are fitted 
by spline interpolation. Each box contains > 1500 s of data. 
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Fig. 3. — Power density spectra (top panels, in units of rms 2 ) with 2-60 keV and the cross 
spectra (middle panels) between 9.8-14.8 keV and 2-4.5 keV on the up turn (box 1), the left 
HB (box 3), the middle HB (box 5), the HB/NB vertex (box 12) and the middle NB (box 
18) from left to right of Interval I. The black circles and the red plus signs represent the 
soft and hard lag, respectively. The corresponding phase lags display in the bottom panels. 
The power density spectra are fitted by a multi-Lorentzian function. See text for the details. 
Each component presents as a dash line and the best fitting results are displayed as grey 
solid lines in top panels, which also appears in Fig. H]-Fig. [5j 
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Fig. 4. — Power density spectra (top panels) and the cross spectra (middle panels) as well 
as the phase lags (bottom panels) on the up turn, the left HB, the middle HB, the right HB 
and the HB/NB vertex of Interval II. The selected energy channels are the same with Fig. 
[3j The PDSs and the cross spectra of the NB and FB are not provided here because the 
data are not dense enough. 
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Fig. 5. — The power density spectra (top panels), the cross spectra (middle panels) and the 
phase lags (bottom panels) on the HB of Interval III. 
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Fig. 6. — QPO centroid frequency, the time lag of QPO and rms of QPO on the HB, as a 
function of S z , of Interval I, from top to bottom. Positive lag means that hard photons lag 
soft photons. The time lag of QPO is averaged between v c — Av/2 and v c + Au/2, where v c 
and Av are the centroid frequency and FWHM of QPO, respectively. 
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Fig. 7. — QPO centroid frequency (top), the time lag of QPO (middle) and rms of QPO 
(bottom) on the HB of Interval II, as a function of S z . 
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Fig. 8. — The QPO centroid frequencies and the time lags of the HB/NB vertex for Interval 
I- IV from left to right, respectively. The HB/NB vertex region is manually selected from the 
relative HID. The time lag is computed between 7.4 - 9.4 keV and 2 - 4.5 keV. 
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Fig. 9. — The NBOs in Interval I as a function of S z . The PDSs on the NB are computed 
from ~ 7 — 60 keV band light curves. The selection regions of NBOs correspond to the box 
14-23 in Fig. El 
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Fig. 10. — Average time lags between 7.4-9.4 keV and 2-4.5 keV in the frequency between 4 
and 8 Hz on the NB of Interval I, where the NBOs around 6 Hz were observed. 
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Fig. 11. — Average time lags on the NB and FB of Interval III. The energy channels and 
frequency range selection are the same with Fig. [TUJ 
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Fig. 12. — The time lag spectrum (left panel) and the rms amplitude spectrum (right panel) 
for the NBOs on the lower NB of Interval I. 
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Fig. 13. — The time lag spectrum (left panel) and the rms amplitude spectrum (right panel) 
for the NBOs on the lower NB of Interval III. 
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Fig. 14. — The time lag spectrum (left panel) and the rms amplitude spectrum (right panel) 
for the NBO on the lower NB of Interval IV. 



